Dark-state suppression and optimization of laser cooling and fluorescence in a trapped 
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We study the formation and destabilization of dark states in a single trapped ^^Sr"*" ion caused by 
the coohng and repumping laser fields required for Doppler cooling and fluorescence detection of the 
ion. By numericaUy solving the time-dependent density matrix equations for the eight-level system 
consisting of the sublevels of the 58^51/2, 5p^Pi/2, and Ad^D^/2 states, we analyze the different 
types of dark states and how to prevent them in order to maximize the scattering rate, which is 
crucial for both the cooling and the detection of the ion. The influence of the laser linewidths 
and ion motion on the scattering rate and the dark resonances is studied. The calculations are 
then compared with experimental results obtained with an endcap ion trap system located at the 
National Research Council of Canada and found to be in good agreement. The results are applicable 
also to other alkaline earth ions and isotopes without hyperflne structure. 

PACS numbers: 32.80.Xx, 42.50.Gy, 32.70.Jz, 37.10.Rs 



I. INTRODUCTION 

A single ion in a radiofrequency (RF) ion trap can 
be laser cooled to a low temperature and confined to 
a region in space with dimensions less than an optical 
wavelength. It thus constitutes a pure and well isolated 
quantum system that lends itself to applications where 
control and isolation from the environment are called for: 
Narrow optical transitions in trapped ions can be used to 
realize optical frequency standards [1, 2], and ion traps 
are also promising systems for the implementation of a 
scalable quantum computer [3, 4]. 

We study even isotopes of the alkaline-earth-metal ions 
that have no nuclear spin and thus no hyperfine struc- 
ture, most notably 24Mg+, ^°Cs+, ^^Sv+ , and ^^^Ba+ , 
but also including less abundant isotopes of these ele- 
ments. Numerical calculations and a comparison with 
experimental results are performed for ^*Sr^, the lowest- 
lying energy levels of which are shown in Fig. 1. The 
'^Si/2— transition is used to Doppler-cool the ion. 

The narrow ~ ^^5/2 electric quadrupole transition 
is used as the clock transition in optical clocks, for sub- 
Doppler cooling using resolved sideband cooling [5], and 
as a qubit candidate for quantum computing [6]. The 
occurrence of a transition to the state is detected 

using Dehmelt's electron shelving technique [7, 8], i.e., as 
a dark period in the fluorescence caused by the cooling 
laser. The ^Pi/2 excited state has a finite probability of 
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decaying to the metastable ^l?3/2 state. This requires a 
repumping laser tuned to the "^0^12 ~ transition 
that returns the ion to the cooling cycle. 

A high scattering rate is crucial for both Doppler cool- 
ing and fluorescence detection of the trapped ion. It is 
compromised by the formation of dark states, which can 
be angular momentum eigenstates or superpositions of 
Zeeman sublevels of the ^S'1/2 and/or ^I?3/2 states (co- 
herent population trapping, CPT [9, 10]). CPT reso- 
nances between the ^6*1/2 and the states can be 
experimentally observed by sweeping the detuning of one 
of the laser beams [11-15]. It is also known that opti- 
cal pumping into dark superposition states between the 
sublevels can occur when the magnetic field is weak 
and that this can be prevented for example by modulat- 
ing the polarization of the repumper [15, 16]. 

Berkeland and Boshier [17] have theoretically studied 
how to destabilize dark states in Zeeman degenerate sys- 
tems using an external magnetic field or by modulating 
the laser polarization. Both methods rely on making the 
dark states evolve in time more rapidly than the optical 
pumping processes. One of the systems they consider 
is the ^S'1/2 — ~ ^1^3/2 system. They recommend 

Rabi frequencies that give a reasonable scattering rate 
without excessive power broadening and point out that 
the repumper should have a positive detuning to tune 
the system away from CPT resonance (the cooling laser 
is red-detuned for Doppler cooling). The magnetic field 
or rate of polarization modulation required to destabilize 
the dark states is then analyzed. 

This paper expands on the results in Ref. [17]. We 
allow the two laser fields to have different polarizations 
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FIG. 1. Lowest energy levels of Sr . Cooling and repump- 
ing transitions are shown as arrows together with their re- 
spective wavelengths. 

and take into account the effect of ttie laser linewidths 
and correlation properties on the CPT resonances. Fur- 
thermore, the effect of the motion of the ion, secular and 
micromotion, is analyzed. Finally, the theoretical results 
are compared to experimental spectra from the ®^Sr^ 
endcap trap at the National Research Council of Canada 
(NRC) [18]. 




m,= -3/2 -1/2 +1/2 +3/2 



FIG. 2. (Color online) The eight-level system. Fine-structure 
state designations are given at the left together with the re- 
spective Lande g factors. Magnetic quantum numbers are 
shown at the bottom. Sublevels are labeled consecutively and 
the allowed transitions between them are shown as arrows to- 
gether with the corresponding relative transition amplitudes 



in Fig. 2. The electric field of the laser field i = c,r is 
written as 



II. THEORY 

The atomic Hamiltonian of the eight-level system 
shown in Fig. 2 is 

2 6 8 

i/a - nwg^\i){i\ + nw,y,^\i){i\ + nuje^\i){il (i) 

i=l 1=3 i=7 

where Wg, Wm, and u,, are the frequencies of the ground 
state \g), metastable state |m), and excited state |e), re- 
spectively. The interaction with an external magnetic 
field B, Vm = — /i.m ■ B, where /Xm is the magnetic dipole 
moment, is included in the atomic Hamiltonian in the 
form of (linearly) Zeeman shifted levels, 

8 

H,^^^^nuj,\i){i\. (2) 
1=1 

Here uJi — uji + gimjj^BB/h, where I is g, m, or e and gi 
is the corresponding Lande factor, mj_,; is the magnetic 
quantum number of the level and fiB is the Bohr mag- 
neton. We have chosen the quantization axis (QA), i.e., 
the z axis, along the magnetic field. The ion is located 
at the origin. It is first assumed to be stationary, later 
its motion is accounted for. 

The system is interacting with two laser fields that 
can have arbitrary time-dependent polarizations. The 
cooling field is slightly red-detuned from the \g) — > |e) 
transition (for Doppler cooling) and induces the transi- 
tions shown as black arrows in Fig. 2, while the repump- 
ing field is near resonant with the |m) — > |e) transition 
and induces the transitions shown as gray (red) arrows 



= y [m^o + /;i(t)u+i + /li(t)u_i] e-»--*-Hc.c., 

(3) 

where £i is the electric field amplitude, uji is the fre- 
quency, c.c. stands for complex conjugate and fq{t) is 
the time-dependent complex amplitude of the polariza- 
tion component q, obtained as the projection of the unit 
polarization vector of the field onto the spherical unit 
vector Uq. These are 

u±i = T^iux ±iuy), uo=U2, (4) 

and describe and tt polarized light, respectively. 

In the electric dipole and rotating wave approxima- 
tions, the interaction Hamiltonians for the two fields be- 
come 

Vcit) = -^[/5W(C7i|7)(l| + C82|8)(2|) 

+ fl,it)Csi\8){l\ 

+ /^i(t)C72|7)(2|]e-''"^*-fH.c. (5) 

and 

Vr{t) = -^[/5W(C74|7)(4| +C85|8)(5|) 

+ /;i(t)(C73|7)(3|+C84|8)(4|) 
+ /iiW(Cr5|7)(5| + C8e|8)(6|)] 
X e-*"'* + H.c., (6) 

where H.c. stands for Hermitian conjugate. The coeffi- 
cients Cji are the relative transition amplitudes of the 
electric dipole transitions and can be obtained using the 
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Wigner-Eckart theorem, which states how a matrix ele- 
ment depends on the mj quantum numbers [19] 



{J'm'j\T^\Jmj) = (-l)-^'-'"'-'(J'||T«||J) 



J' K J 

-m'j q mj 



rpK 

— ^J'm'j,Jmj->- Ji J- 



(7) 



Here is a set of tensor operators of rank k and Tj, j = 

{J'\\T'^\\J) is the reduced matrix element that docs not 
depend on mj, m'j, or q. Using Eq. (7) for the electric 
dipole operator (k = 1; g = for tt and ±1 for 
transitions), one obtains: C^2 = —C71 = 6~^^'^, C72 = 
— Cgi — 3^^/^, C73 — Cs6 = 2~^, C74 = Cs5 = — 6~-^/^, 
and C75 = Cs4 = 2-i3-i/2_ 

The "two-level" Rabi frequencies corresponding to the 
two laser fields are flc = Mf'i/2,5'i/2'^c/^ fir = 
MPi/2.-D3/2'^r/fi, where IJ.l'^,,Lj = {L'j,\\ - er\\Lj) is the 
reduced dipole moment. The Rabi frequencies can be 
assumed to be real. The true Rabi frequency for a cer- 
tain transition |z) — >■ \j) is equal to Cji^lk (k = c, r 
depending on which laser field drives the transition). 
Note that Ref. [17] uses rms Rabi frequencies defined as 



^'rms A-in 



mjrrij 



I p. Prom the sum rule [19] 



mj,mj 



(8) 



we obtain the relation between the rms and the two-level 
Rabi frequencies: firms = S~^^^fl. 

The Liouville-von Neumann equation for the density 
matrix of the eight-level system is 

do 1 do"'^'''^ 

^ = ^[//a,. + KW + mp] + ^. (9) 

Terms oscillating at optical frequencies are eliminated by 
the substitutions 



Pki = Pkie 



-iojct 



l<k<2, 7<l< 



(10a) 

Pkl=Pkie~'"'', 3<k<6,7<l<8, (10b) 
Pki = Pkie-'^"'-"''>', 1 < A; < 2, 3 < / < 6. (10c) 



dt 



herences in the spontaneous emission process [20] 

(2 Je + 1) ^ Cj Cj 
q=-l 

^ Pmc=mg+q,Tn'^=Tn'g+qt 

(12) 

where the transition amplitudes Cj^me,Jgmg are as de- 
fined above. This adds the following relaxation terms to 
the evolution equations 

jj„relax 

^'•'+1 = 2CnCs,^+lAmrp7s, i = 3, 4, 5, (13a) 



dt 

(Jpfelax 

— — — = 2CnCs2Agrprs- 



(13b) 



The optical coherences decay at the rate r/2 (wc; as- 
sume the laser linewidths to be <^ T) and the metastable 
state populations and coherences decay at the rate jm ^ 
r. The transition probabilities for the quadrupolc tran- 
sitions can be obtained from Eq. 7 with k = 2 and the 
branching ratios for decay from the metastable state are 

then Amg = \Cmg\VEg\Cmg\'': ^62 = Asi = 1/5, 

-461 = A32 = 4/5, A52 = A41 = 2/5, and A51 = A42 = 
3/5. 

Wc have assumed that there are no external processes, 
such as background gas collisions, that cause transitions 
between the ground or metastable m j levels or dephasing 
of the ground and metastable coherences. The dephas- 
ing rate of the coherences between the ground and the 
metastable levels depends on the laser linewidths [10] 



'Jg.rn 



2AWc,r) 



(14) 



where Aujc and Aoj^ are the FWHM linewidths of the 
cooling and rcpumper lasers and Awc.r is the cross- 
correlated linewidth. For perfectly cross-correlated lasers 
the linewidths cancel, whereas for completely uncorre- 
lated lasers AlOc,t = and the linewidths add up. In the 
low-intensity limit, the linewidth of the CPT resonances 
is 27g,„. 



A. RelcLxation 

The relaxation terms dp'^^^'^^/dt in Eq. (9) are added 
phenomenologically. The excited state populations decay 
at the rate F. The excited state decays to the ground and 
metastable states with the probabilities Ag and Am = 
1 — Ag and the decay probability for a certain magnetic 
sublevel transition, e.g., |8) to |6) can then be calculated 
as 

\CS6? + \CS5\' + \CS4\^' ^ ' 

The excited state coherence pis decays at the rate F 
and can be partially transferred to the ground state co- 



B. Solution of density matrix equations 

The time-dependent density matrix is solved by direct 
numerical integration of Eq. (9) with the initial condi- 
tion pii = P22 = 1/2, i.e., a thermally populated ground 
state. This describes the situation well when the cooling 
and repumping beams are turned on after a probe-laser 
interrogation on the ^5*1/2 — ^£'5/2 quadrupole transi- 
tion. When the ion is first trapped, the situation is more 
complicated, as it is hot and there might be another far- 
detimcd cooling field present [21]. However, also in this 
case the treatment should be valid for the "final" state, 
i.e., when the ion has cooled down and reached quasi- 
steady state (when the only time dependence is due to 
the possible polarization modulation of the repumper). 
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TABLE I. Numerical parameters for ^^Sr^ 



Parameter 


\aiuo 


llrf. No. 




3.078 eao 


[22] 


A*^'l/2.03/2 


3.112 eao 


[22] 


r 


135.58 X 10*^ s"^ 


[22] 


A, 


0.9444 


[22] 


7rr> 


2.30 s"^ 


[26] 



As the density matrix describes ensemble averages, the 
time-dependent solution should be interpreted as the av- 
erage over an infinite number of experimental cycles for 
a single trapped ion. 

If the polarization amplitudes /* are time indepen- 
dent, the steady-state solution can be solved by setting 
the time derivative to in Eq. (9). If they are time- 
dependent, we must integrate Eq. (9) until quasi-steady 
state is reached and average the density matrix over one 
polarization modulation period, which is considerably 
more computer time consuming. 

The (quasi-)steady-state value of the scattering rate 
is the most relevant parameter for comparison with ex- 
periments. To collect the fluorescence of the single ion, 
it is customary to use a bandpass filter at the cooling 
wavelength to minimize problems with background light 
and scattering. We therefore define the scattering rate 
as = AgV^p'jT + pss). 



III. EXPERIMENTAL PARAMETERS FOR 

88SJ.+ 

Table I summarizes the numerical parameters used for 

the ^*Sr+ ion in this paper. We use the dipole matrix 
elements, spontaneous decay rate, and branching ratio 
calculated by Jiang et al. [22]. Their values are in good 
agreement with the experimental lifetime in [23] 

and with the calculated lifetime and oscillator strengths 
in [24] . Gallagher's [25] experimental oscillator strengths 
and branching ratio are often cited. However, his oscilla- 
tor strength for the — ^-Pi/2 transition agrees with 
the results in [22] , but he acknowledges that his value for 
the ^£'3/2 — transition probably is too large, which 
is why we chose to use the values of Jiang et al. [22]. 

Using li = eoc£f/2, we can relate the two- level Rabi 
frequencies to the laser intensities 



(15) 



where is shorthand notation for the corresponding 
reduced dipole moment. This gives the numerical re- 
lations ric/r = (/c/39.8 mWcm-2)i/2 and n./T = 
(/r/39.0 mWcm~^)^/^. The laser beams can be carefully 
adjusted, so we assume that the ion experiences the peak 
intensity of the Gaussian beam, Ij = 2Pi/TTW^, where Pi 
is the power and Wi is the waist. 



Berkeland and Boshier [17] recommend using the two- 
level Rabi frequencies Oc « Oj. « 3~^/^r in order to ob- 
tain a reasonable scattering rate with little power broad- 
ening. However, considerably higher cooling laser Rabi 
frequencies have been used in experiments at both the 
NRG (r^c ~ 20r) [27] and the National Physical Labora- 
tory (Qc « 7r) [28]. 

IV. TYPES OF DARK STATES 
A. Ground-level dark states 

Berkeland and Boshier [17] state that for a transition 
between states with half-integer angular momentum and 
Jg = Je, such as the cooling transition considered here, 
there is one dark state for circular polarization only. Let 
us elaborate on this. 

Figure 3 shows the scattering rate as a function of time 
for different cooling laser polarizations. For a circularly 
polarized cooling beam in zero magnetic field, we can 
choose the QA such that the polarization is a~^. This 
polarization only drives the |1) — > |8) transition and 
the scattering rate decays rapidly as the population is 
pumped into the dark state |2) at the optical pumping 
rate Fop = 0.045F [solid (black) curve in Fig. 3]. We have 
earlier derived an expression for the optical pumping rate 
in a three-level system with one dark lower level [29] 



■ op 



(16) 



Here An is the spontaneous decay probability to the non- 
coupled state. If we neglect the metastable state, Eq. (16) 
can be applied to the effective three-level system consist- 
ing of |1), |2), and |8) (A„ = As2, C = Cgi), giving 
Fop = 0.04ir in reasonable agreement with the numeri- 
cal result above. If we apply a magnetic field in the di- 
rection of the cooling beam, 1 2) is still a dark eigenstate 
and the Zeeman shift only changes the detuning. 

If we choose the QA in another direction, e.g., in the 
X direction when the beam is in the z direction, the po- 
larization in the QA frame is — 2^-'^/^[uo — z2^^/^(u+i — 
u_i)], i.e., a combination of tt and a (cr refers to light 
that has a linear polarization orthogonal to the QA and 
that thus induces (T+ and cr" transitions). The dark state 
is now the superposition state 



lDSi,2) = ^(|l)-i|2)), 



(17) 



for which (7|Vc|DSi,2) = (8|K|DSi,2) = 0. If we apply a 
magnetic field in the direction of the QA (u,.), |DSi_2) is 
no longer an energy eigenstate, but the dark state evolves 
in time 



DSi,2)(t) = ^^(e'9s''BBt/2fi|^^ _ ■g-ig,^BBt/2fi|2))_ 



V2 



(18) 
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FIG. 3. (Color online) Time-dependent scattering rates for 
different cooling polarizations. A polarized cooling beam 
rapidly pumps the population into the dark state |2) [B = 0; 
solid (black) curve]. In a transversal magnetic field, the dark 
state is a superposition state, which is destabilized by the field 
[B = 350 ^Tux ; dashed (red) curve] . A combination of tt and 
cr"*" polarization creates a partially dark state [B = 0; dash- 
dotted (green) curve]. For a linearly polarized cooling beam, 
there are no ground-level dark states [B = 0; dotted (blue) 
curve]. For all curves, Qc = = F, 5c = -F/2, = F/2, 
79, m = F/fO, and repumper polarized along Ui,. 



When the field is strong enough that the evolution rate of 
the dark state gg^^B/h is higher than the optical pump- 
ing rate Fop (corresponding to « 35 /iT for the param- 
eters in Fig. 3), the system cannot follow and the dark 
state is destabilized. Another way to explain this is that 
the optical pumping rate Fop determines the width of the 
ground states, and the Zeeman shift of the ground lev- 
els tune the A system out of Raman (CPT) resonance. 
When ggfisB/h « lOFop {B « 350 /uT), the dark state 
is completely destabilized; see the dashed (red) curve in 
Fig. 3. If the field is increased further so that the Zeeman 
shifts approach F, the cooling laser is tuned away from 
resonance and the scattering rate decreases. 

If the polarization is a combination of tt and (this 
can be achieved by an elliptically polarized beam travel- 
ing at an angle to the QA), the A system |1) o |8) O |2) 
is formed, but it is not closed, as the |1) -;-> |7) transition 
is also driven. Hence the contrast of the corresponding 
CPT resonance is less than unity, i.e., the scattering rate 
decays rapidly, but not to [dash-dotted (green) curve 
in Fig. 3]. 

If the cooling laser is linearly polarized (tt, a, or com- 
bination thereof in the QA frame), there is no dark state 
among the ground levels. The slow decay of the dotted 
(blue) curve in Fig. 3 is due to dark states among the 
metastable sublevels, which are considered in Sec. IV C. 
In the following, we therefore only consider a linearly 
polarized cooling beam. 



B. Ground- metastable- level dark states 

The ground-metastable-level dark states are CPT su- 
perposition states of \g) and |m) sublevels. As the two 
legs of the A systems are driven by separate lasers, the 
CPT resonances can be detected by tuning one of the 
lasers over the Raman resonance [11-15]. 

The width and contrast of the ground-metastable CPT 
resonances depend strongly on the ground-metastable co- 
herence dephasing rate 7g,m, Eq- (14). Siemers et al. [13] 
used narrow-band lasers ("a few kHz") and deduced a 
coherence dephasing rate of 27r x 12 kHz, in agreement 
with the assumption of noncorrelated lasers (the actual 
CPT resonances were wider due to power broadening). 

The dephasing due to the laser linewidths was not ex- 
amined in the calculations in Ref. [17]. Figure 4 shows 
ground-metastable CPT resonances for parameters cor- 
responding to curve A in Fig. 7 of Ref. [17], but for dif- 
ferent dephasing rates 7g,m- As this rate approaches F, 
the CPT resonances disappear completely. The magnetic 
field of 150 /xT is large enough to separate the four res- 
onances as well as to destabilize metastable-level dark 
states (see Sec. IV C). The CPT resonances have finite 
linewidths even for ^g^m = due to power broadening, 
which is illustrated in Fig. 5. Note that the blue side of 
the ^5*1/2 — ^ Pi/2 transition cannot be scanned in a real 
ion trap experiment as the ion is lost due to heating. 




0.2 0.4 0.6 0.8 1 
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FIG. 4. (Color online) CPT resonances for different values 
of 7g,m/r: [solid (black) curve], 0.01 [dashed (red) curve], 
and 0.1 [dotted (blue) curve]. Q.c = Vlr = V /2, 5, = r/2, 
^BB/h = r/10 (B ^ 150 pT). 



As the cooling laser is red-detuned (typically 5c ~ 
~F/2), ground-metastable CPT resonances can be 
avoided by tuning the repumper slightly to the blue. If 
the repumper is modulated so that sidebands are cre- 
ated, care must be taken that no sideband is at Raman 
resonance either, as discussed in Sec. VI. 
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FIG. 5. (Color online) Power broadening of the CPT res- 
onances. Normalized (at Sc = F) scattering rate spectra 
for different Rabi frequencies fie ~ ^t'- O.IF [solid (black) 
curve], 0.3r [dashed (red) curve], and F [dotted (blue) curvee]. 
79,m = 0, = F/2, fisB/h = F/IO (S « 150 fiT). 



FIG. 6. (Color online) For B = 0, the scattering rate is 
independent of the repumper polarization [tt, a, or ct"*"; solid 
(black) curve]. For a a polarized repumper, a field of 100 fiT 
{2gml-iBB /h ~ 0.09r) destabilizes the dark states completely 
[dashed (red) curve], fic = fir = F, Sc = -F/2, Sr = F/2, 
79, m = F/lO, and cooling laser a polarized. 



C. Metastable- level dark states 

The \Jm = 3/2) — ;> |Je = 1/2) repumper transition 
has dark states for all laser polarizations. In zero field, 
the scattering rate is practically independent of the re- 
pumper polarization [see the solid (black) curve in Fig. 6]. 
The optical pumping rate is Fop ~ 0.002r and again 
Eq. (16) gives a good estimate if we assume both fields 
to be linearly polarized so that we can use C = Cs2 
and An ~ A^q. This optical pumping rate is lower than 
the optical pumping rate into ground-level dark states 
(Sec. IV A) by roughly a factor of A„i/Ag as one would 
expect. The actual dark states depend on the repumper 
polarization. For tt polarization, |3) and |6) are dark 
states, for ct^, |5) and |6) are dark, and for a polariza- 
tion, the dark states are superpositions of |3) and |5) and 
of |4) and |6). 

When the repumper is a polarized, the dark states can 
again be destabilized using a magnetic field. The relation 
between the optical pumping rate and the required field 
is not as simple as in Sec. IV A because the decay prob- 
ability A„i is low. However, a field of 100 /iT (dark state 
evolution rate 2g,nfJ'BB /h ~ 0.09F) is enough to fully 
destabilize the dark states for the parameters in Fig. 6 
[dashed (red) curve]. 



V. POLARIZATION ANGLES 

Let us summarize the requirements on the laser polar- 
izations. In Sec. IV A we showed that the cooling laser 
beam should be linearly polarized. If we are far from 
Raman resonance, its function, from the point of view 
of dark states, is merely to pump population into the 
metastable levels, where the repumper can then build 



up coherences (dark states) unless these are destabi- 
lized. Hence the effect of the angle between the mag- 
netic field and the linear cooling laser polarization is 
negligible, typically only a few percent for the param- 
eters considered here. On the other hand, near Raman 
resonance the number, position, and amplitude of the 
ground-metastable CPT resonances depend strongly on 
the polarization of both laser fields, but this is a regime 
we try to avoid. 

The polarization angle of the repumper is crucial for 
the magnetic-field destabilization to work (see Fig. 7) . It 
works only when the dark states are superpositions, i.e., 
when the repumper contains both cr"*" and a~ . For all 
other polarizations (tt, cr+, , combination of tt and ct"*" 
or (7^, but not both) there is at least one dark eigenstate 
into which optical pumping occurs. 



VI. POLARIZATION MODULATION 

Dark state destabilization using an external magnetic 
field has been discussed in Sec. IV. However, there are 
applications where a strong magnetic field is not desired, 
in particular ion clocks, where typically fields of only a 
few microteslas are used. 

A simple and commonly used way to modulate the 
repumper polarization is to pass the beam through an 
electro-optical modulator (EOM) at a 45° angle to the 
optical axis. The resulting polarization (with beam di- 
rection as QA) is 



UEOM(t) = (1 + 



V2 



-iVEOM(t) 



-iVEOM(t), 



)u+i + (1 - ie' 



-i¥'EOM(t) 



)u_i 
(19) 
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FIG. 7. (Color online) Scattering rate as a function of the 
angle between the magnetic field and the linear repumper po- 
larization when dark states are destabilized using B = 50 /iT 
[solid (black) curve] and B = 100 /iT [dashed (red) curve], 
fie = f^r = r, Sc = -r/2, 5r = r/2, 7g,™ = r/lO, and cooling 
laser a polarized. 



where the phase retardation is [17] 

<PEOM{t) = ^$EOm(1 - COSWEOM^)- (20) 

If the modulation frequency weom = or the modula- 
tion amplitude $eom = 0, the resulting polarization is 
the stationary linear polarization incident on the EOM, 
(Uj. + Uy)/\/2. In the following, we consider only the 
EOM modulation technique. Modulation techniques us- 
ing acousto-optic modulators (AOMs) give similar results 
[17]. 

Figure 8 shows that the dark state destabilization due 
to polarization modulation is very similar to magnetic 
field destabilization when the dark state evolution rate is 
the same (cf. Fig. 6), except for some oscillations at the 
modulation frequency weom- The two effects also sum 
up almost fully [see the dotted (blue) curve in Fig. 8]. 

Figure 9(a) shows the quasi-steady-state scattering 
rate as a function of the modulation frequency weom- 
As shown in Fig. 8, the dark-state destabilization is ef- 
ficient already at (^eom ~ O.IF for the parameters used 
here. As expected, there are dips at F and r/2, where 
the first- and second-order red sidebands of the repumper 
are at Raman resonance. Assuming we have a fixed EOM 
modulation frequency, this must be taken into account 
when choosing the repumper detuning. The dips be- 
come more pronounced if the Rabi frequencies are in- 
creased or if 7g^„i is decreased. Figure 9(b) shows that 
for cjeom = r/4, the scattering rate is maximized for the 
modulation amplitude $eom ~ lAn. At lower ampli- 
tudes the polarization variation is not sufficient to fully 
destabilize the dark states, and at higher amplitudes an 
increasing part of the repumper intensity is transferred 
into sidebands further from resonance. 
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FIG. 8. (Color online) Scattering rate for stationary a re- 
pumper [solid (black) curve] and polarization modulated ac- 
cording to Eqs. (19) and (20) at wbom = 27r x 2 MHz f» O.OQF 
[dashed (red) curve]. Modulation at 1 MHz plus a 50-/iT field 
along Ux gives approximately the same result [dotted (blue) 
curve]. $EOM = TT, fie = a = r, Sc = -r/2, Sr = F/2, 
7s, m = r/10, and cooling laser u^; polarized. 
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FIG. 9. Quasi-steady-state scattering rate (a) as a function 
of modulation frequency weom for $eom ~ tt and (b) as a 
function of modulation amplitude "3?eom for ljeom = r/4. 
Qc = ^ F, Sc = -F/2, 5, = F/2, 79,™^ = F/10, coohng 
laser u^, polarized, and repumper as in Eqs. (19) and (20). 



VII. ION MOTION 

The ion undergoes both secular (thermal) motion at 
the trap frequencies and micromotion at the applied trap 
drive frequency. If the ion is not significantly perturbed 
from the RF voltage symmetry node of the trap by stray 
fields so that the micromotion amplitude is kept low, the 
Doppler cooling is not compromised and the ion experi- 
ences an approximately constant laser intensity near the 
beam focus. Under these circumstances, the motion can 
be accounted for by considering the rest frame of the ion 
and including the oscillations in the phases of the laser 
fields. 

The equations of motion of a single ion in a quadrupole 
trap are given by the Mathieu equation and the first- 
order solution in the trap coordinate system XYZ = 
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RxRyRz is [30] 

Ri{t) « [Roi + i?soc,i cos (Wsec.ji + 6*800,1)] + ^ '^^'^ ' 

(21) 

Here Rscd is the secular motion amplitude, which can be 
obtained from the temperature of the ion, i?Ki = k^Ti 
\rnRl^^^ujl^^^i [30], uJscc,i is the trap frequency, and 6*860, i 
is the phase, determined by the initial position and mo- 
mentum of the ion. The trap drive frequency is $7 and qi 
is a parameter that depends on the trap geometry, drive 
voltage, and ion mass. i?oi is the displacement of the av- 
erage position of the ion that is caused by static electric 
fields and that results in excess micromotion. 

We consider a geometry, where the endcap trap Z axis 
is at an angle a = 35.3° to the horizontal plane. The 
possible beam directions are vertical (z) and horizontal 
at ±45° to the horizontal projection of the trap axis [x 
and y). In the beam coordinate system, the ion motion 
is then given by 

r^{t) = [Rx{t) - smaRvit) - cos a Rz{t)] (22a) 
v2 

ry{t) = [Rx{t) + sin aRvit) + cos aRz{t)]{22h) 
v2 

(t) = - cos a Ry (t) + sin a Rz (t) . (22c) 

This geometry is chosen to make the projection of each 
beam direction on the trap Z axis equal in magnitude 
(cosa/-\/2 — sin a). In the rest frame of the ion, the 
phase of a laser beam traveling along the u; axis is then 
given by exp {-i[ujjt - kjri{t) - (f>ji]} {j = c, r). The 
phase is not relevant. 

Figure 10 illustrates the effect of ion motion using nu- 
merical parameters typical for an endcap trap: qz = 
-2qx,Y = 0.4, u}sec,z = 2uJsec,x,Y = 2TT X 2 MHz, and 
17 = 27r X 12 MHz. The solid (black) curve is the scatter- 
ing rate for a stationary ion, T = 0, equal to the dashed 
(red) curve in Fig. 8 (note the different horizontal axis), 
and shows oscillations at the frequency weom only. At fi- 
nite temperatures, there are additional oscillations at the 
trap frequencies Ll)scc.x,y and uJscc,z and at the drive fre- 
quency S7, and the time dependence of the scattering rate 
is rather complex. The oscillation amplitudes increase 
with increasing temperature; see the dashed (red) and 
dash-dotted (green) curves in Fig. 10, which correspond 
to temperatures of 0.5 niK (Doppler limit for ^^Sr"*") and 
5 mK, respectively. However, as the density matrix for- 
malism describes an ensemble average over the internal 
degrees of freedom, it is not physically meaningful to se- 
lect fixed values for the external degrees of freedom, i.e., 
the secular motion phases 9i that depend on the initial 
position and velocity of the ion. Instead, we must aver- 
age over a large number of solutions with random phases 
9i. The dotted (blue) curve in Fig. 10 shows that an av- 
erage over 250 solutions with T — h mK is very close to 
the curve for a stationary ion. 

Intuitively, one expects the ion motion to have the 
most significant effect in the vicinity of a narrow CPT 



resonance. The left-hand inset in Fig. 10 shows CPT 
resonances for 7g,m = 0. For a stationary ion (solid 
black curve), the scattering rate vanishes at CPT res- 
onance, whereas the resonance is significantly "rounded 
off" in the T — 5 mK case (dashed red curve). However, 
the right-hand inset in Fig. 10 shows the corresponding 
curves for 7g,m = 0.04r (the value that is used in the 
comparison with experiments in Sec. VIII) and now the 
difference between the curves is minute. For the inset 
curves, the ensemble average was evaluated as a time 
average once quasi-steady-state was reached, which was 
particularly convenient as the involved frequencies were 
chosen to be multiples of u1soc,x,y- 

Thus, unless one is studying CPT resonances with a 
very low coherence dephasing rate "fg,rm one can closely 
approximate the scattering rate without accounting for 
ion motion in cases where the ion motion is of a low 
amplitude and micromotion effects do not dominate. 
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FIG. 10. (Color online) Scattering rate for a stationary ion 
[T — 0; solid (black) curve] and at temperatures 0.5 mK 
[dashed (red) curve] and 5 mK [dashed-dotted (green) curve]. 
The average over 250 curves with T = 5 mK and random 
secular motion phases 6i approaches the stationary case [dot- 
ted (blue) curve]. Left inset: CPT resonances for a stationary 
ion [solid (black) curve] and at T = 5 mK [dashed (red) curve] 
for 7g,m ~ 0. Right inset: Corresponding CPT resonances for 
7a, m = 0.04r. Parameters not mentioned are the same as for 
the dashed (red) curve in Fig. 8. 



VIII. COMPARISON WITH EXPERIMENTS 

A. Endcap trap setup 

We now compare the theoretical results to experimen- 
tal spectra measured using the ^^Sr"*" endcap trap at NRC 
[18]. Briefiy, the experiment employs a single ion of 
**Sr^ held in an endcap trap [31, 32]. The trap char- 
acteristic dimensions have endcap electrodes of radius 
rg = 0.25 mm with separation 2zo = 0.541 ± 0.015 mm 
made of molybdenum. The shield electrodes were kept 
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at ground potential, had a diameter of 2r2 = 2 mm, and 
were made from tantalum tubing. The entire trap striic- 
ture was composed of nonmagnetic materials and housed 
in an ultra- high- vacuum environment, where the back- 
ground pressure was 10 nPa, composed primarily of H2. 
At these pressures, the mean collision times are of the or- 
der of 10^ s [33] and do not contribute to the dephasing 
rates of the observed line shapes. The ion trap is operated 
at a trap frequency of = 27r x 14.4 MHz with a voltage 
amplitude of Vq = 200 V. The trap secular frequencies 
were Wsec,z/27r = 2.3 MHz and a;sec,x,r/27r = 1.2 MHz in 
the radial direction, with the difference between the two 
radial frequencies being about 20 kHz. The ion is loaded 
into the trap using a weak effusive source of Sr atoms and 
photoionized using 461- and 404- nm laser beams which 
are turned off following the loading process [21]. In this 
way, micromotion due to patch potentials of deposited Sr 
is reduced to a minimum and only small adjustments of 
the trap compensation electrodes are necessary between 
loadings. The trap is enclosed in a double-layer mag- 
netic shield and the residual field at the trap center is 
B = (-1.22u^ - 0.70uy 1.47u^) yuT (|B| = 2.03 /zT). 

In the experiment, both beams propagate in the Uj, 
direction. The cooling laser is linearly polarized in the 
Mz direction, whereas the repumper polarization is mod- 
ulated like = 2-i/2(u,; + e-*'^EOM(0u^). The EOM 
frequency is 12 MHz (weom = 0.556r) and the ampli- 
tude is $EOM = (1-1 ± 0.2)77. Typical beam parameters 
at trap center are = 16 ± 1 /xm and = 50.5 ± 2 /im. 
The linewidth of the repumper (a diode pumped fiber 
laser) is specified to below 10 kHz and can be neglected. 
The repump laser is stabilized in its drift by stabilization 
to a transfer optical cavity whose length is controlled 
using a 633-nm polarization-stabilized HeNe laser sys- 
tem which maintains the repumper laser within 2 MHz 
over extended periods of time. The cooling laser system 
has been described in a previous publication [34]. The 
laser is a 422-nm diode laser system which is stabilized 
in its short-term frequency fluctuations to the side of a 
transmission fringe of an optical cavity and then con- 
trolled in its long-term drift by locking to the saturated 
absorption of the hs'^Siji — 6p^Pi/2 fine in ^'Rb. The 
linewidth of the cooling laser has been determined by 
two different independent techniques to be 2.4±0.8 MHz 
[34] and 0.6 ± 0.3 MHz. A ground-metastable coherence 
decay rate 'jg^m = 0.04r, corresponding to a linewidth 
of 1.7 MHz, is used in the current comparison, as it 
gives the best agreement between experiments and the- 
ory. The cooling laser is frequency referenced to the ^^Rb 
5s^Sy2 {F = 2)- 6p^Pi/2 {F = 3) transition [34, 35], 
which provides the absolute frequency scale for the cool- 
ing laser line shapes. In the repumper line shapes, the 
absolute frequency scale is determined from the position 
of the CPT resonances. 

The laser Rabi frequencies are determined by measur- 
ing the beam powers and waists, giving an uncertainty 

of approximately 20-25%. However, the slightest mis- 
alignment of the beams affects the Rabi frequencies ex- 



perienced by the ion, so the actual uncertainty is larger. 
The ion trap fluorescence is optimized at low laser in- 
tensities with neutral density filters in the beam paths. 
The filters are then removed for the measurements at 
high Rabi frequencies and it is possible that this alters 
the beam alignment so that the true Rabi frequencies are 
lower than estimated in these cases. 

Ion fluorescence is detected by an / = 27 mm (/- 
number //# = 0.9) aspheric collection lens located in 
the vacuum chamber which focuses the light onto a pin- 
hole in front of a photon counting photomultiplier sys- 
tem. Detected count rates at optimum fluorescence are 
typically at 10 000 cps. The detected fluorescence rate 
is sensitive to scattered cooling laser light, mainly from 
the trap electrodes. This background has not been re- 
moved from the experimental data. However, the level 
of parasitic light in the current experiment is typically 
less than 1/50 of the observed single ion fluorescence un- 
der optimum conditions. The relation between the de- 
tected fluorescence rate and the total scattering rate is 
not known and the experimental data are hence given in 
arbitrary units. A second viewport opposite the photon 
coimting photomultiplier has a commercial camera lens 
system (//# = 4) which images the trap central region 
on a photon counting CCD camera system. This allows 
the initial optimization of beam alignment and ion flu- 
orescence to be performed. In addition, it is used as a 
first stage of minimizing the micromotion using external 
compensation electrodes in the radial plane and axial dc 
bias of the endcap electrodes in the axial direction. Final 
optimization of the ion micromotion is achieved by the 
measurement and optimization of the observed sideband 
features in the spectra of the ion probed on the reference 
ion clock transition at 445 THz. In this way the total 
Stark and time dilation perturbations of the single ion 
are maintained below the 2 x lO"^'^ fractional shift for the 
reference transition at 445 THz and micromotion effects 
are deemed negligible for the current line-shape scans. In 
addition, scans of the ion secular motion sidebands have 
placed an operating temperature of 2 ± 1 mK for all trap 
canonical directions under optimum laser cooling, thus 
well satisfying the condition for Doppler free spectra of 
the ion line shapes. Departures from the ideal low ion ki- 
netic temperatures are anticipated when the laser powers 
cause low fluorescence scattering and detunings are near 
the ion line center. 



B. Results 

Figure 11 shows the scattering rate as a function of the 
EOM modulation amplitude as in Fig. 9(b). The only ad- 
justed parameter in the comparison was the ratio of the 
observed detected count rate to the calculated absolute 
scattering rate, which was set at 8150. Other parameters 
were based on the measured experimental conditions. As 
shown, the observed dependence of fluorescence matches 
the calculated relation very well confirming optimum fiu- 
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FIG. 11. (Color online) Measured [filled (black) circles] and 
calculated [solid (red) curve] scattering rate as a function 
of EOM modulation amplitude. Qc/T = 1.7, fir/r = 0.77, 

Sc/r = -0.74, s,/r = 0.7. 

Figure 12 shows line shapes obtained by scanning the 
cooling laser detuning for increasing Rabi frequencies. 
The repumper was red-detuned in order to make the CPT 
resonances visible. Only the red half of the spectrum 
can be measured and in some of the experimental curves 
one can see how the scattering rate drops close to Sc = 
as the laser cooling is compromised. The agreement 
between the experimental and the calculated line shapes 
is very good, except for the highest Rabi frequency line 
shape. In this case, both the optical line and the CPT 
resonance are considerably more power broadened in the 
calculated curve, indicating that the experimental Rabi 
frequencies are lower than estimated. This is probably 
due to beam misalignment as mentioned above. In the 
second lowest [dashed (red)] curve, the CPT resonance 
is offset, as the repumper detuning has drifted slightly. 
Figure 13 shows the corresponding line shapes for more 
normal operating conditions with a repumper detuning of 
O.OQF. Again, lower fluorescence levels are observed near 
line center, presumably due to ion heating. If a narrower 
linewidth is desired, even lower Rabi frequencies have to 
be used at the sacrifice of the ion scattering rate. 

Figure 14 shows line shapes obtained by scanning 
the repumper frequency for two cooling laser detunings. 
Again, the agreement is very good considering the un- 
certainties of the experimental parameters. The multiple 
coherence dips in the spectra arise from the different or- 
ders of the EOM frequency falling in Raman (CPT) reso- 
nance. Again it should be noted that a significant change 
in the observed line shape is seen for small changes in 
cooling detuning which is well matched by the calcula- 
tions. 

Figure 15 shows repumper line shapes for increasing 
cooling laser Rabi frequencies. The two-peak spectrum 
that occurs for the highest flc values can be explained us- 
ing dressed states [36] . The dressed states are the eigen- 
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FIG. 12. (Color online) (a) Measured and (b) calculated scat- 
tering rate as a function of cooling laser detuning for dif- 
ferent Rabi frequency pairs [fic/F, Slr/F]. From bottom to 
top: [0.36, 0.39] [solid (black) curve], [0.64, 0.70] [dashed (red) 
curve], [1.13, 1.25] [dash-dotted (green) curve], [1.60, 2.22] 
[dotted (blue) curve], and [3.58, 3.22] [solid (magenta) curve]. 
Sr/r = -0.53. 



states of the interacting ion {\g) and |e) levels) plus cool- 
ing laser system when the laser is strong {Cs2^c ^ P)- 
The left-hand side of Fig. 16 shows the uncoupled states 
of the ion-photon system. The energy level scheme con- 
sists of a ladder of manifolds, separated by the cooling 
laser frequency Wc, two of which are shown in the fig- 
ure. Each manifold consists of two closely spaced {\Sc\) 
states \g,N) and |e,7V — 1), where the first letter refers 
to the state of the ion and the second to the number 
of cooling laser photons. When the ion-laser interaction 
is taken into account, the strong cooling laser couples 
the ground and excited states and the new eigenstates 
are the two dressed states |l(iV)) and |2(7V)) shown at 
the right in Fig. 16. These are both superpositions of 
\g) and |e) and are separated by the generalized Rabi 
frequency fl' = a/ C^2^c + ^c- The frequencies of the 
two dressed states in relation to the unperturbed excited 
state are Aa;i.2 — ((5c±ri')/2. Since the \g) — |e) system 
is not closed, the population will be optically pumped 
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FIG. 13. (a) Measured and (b) calculated scattering rate as 
a function of cooling laser detuning for S-^/T = 0.09. Q,c/V = 
1.60, a/r = 2.22. 
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FIG. 14. (Color online) (a) Measured and (b) calculated scat- 
tering rate as a function of repumper detuning for two cool- 
ing detunings: S^/T — —1.30 [solid (black) curve] and 5c/T — 
-0.74 [dashed (red) curve]. flc/T = 1.69, fi^/r = 0.72. 



into \m) except when the repumper is resonant with one 
of the dressed states. For the highest ilc hne shape in 
Fig. 15, the dressed state frequencies are Awi = l.OF 
and Auj2 — — 1.9F. Thus the simple dressed state model 
predicts the peak positions very accurately. In addition, 
each dressed state peak has a CPT resonance at the very 
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FIG. 15. (Color online) (a) Measured and (b) calculated scat- 
tering rate as a function of repumper detuning for increasing 
cooling Rabi frequencies f^c/F: 0.53 [solid (black) curve], 1.19 
[dashed (red) curve], 3.00 [dash-dotted (green) curve], and 
6.71 [dotted (blue) curve]. 4/r = -0.92, STir/F = 0.72. 



center. In practice, such high applied cooling intensities 
should not be used if low ion kinetic temperatures are 
desired, as the high intensities induce power broadening 
and a reduction in the obtainable minimum laser cooling 
temperatures [37]. 

If instead the repumper has a relatively high Rabi fre- 
quency [Fig. 17(a)], a dispersive line shape is obtained, as 
confirmed in the experimental spectrum. The "bright lat- 
eral part" next to the dark resonance was first observed 
by Alzetta et al. [38] and occurs when the two Rabi fre- 
quencies are considerably different and the detunings are 
non-zero. It is caused by stimulated Raman scattering 
between the ground and the metastable states [39] and 
has been studied also in single trapped ions [11, 13]. 

Figure 17(b) shows a cooling laser line shape for typ- 
ical trap operation parameters. The repumper is blue- 
detuned, but the exact detuning is not known. Strong 
ion fluorescence is obtained with the linewidth approach- 
ing the natural width value. The FWHM of the line is 
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FIG. 16. (Color online) Dressed state model: uncoupled levels 
at the left and dressed states at the right. 




FIG. 17. (Color online) (a) Measured [solid (black) curve] and 
calculated [dashed (red) curve] scattering rate as a function 
of repumper detuning for a high repumper Rabi frequency 
a/r = 5.90. 5c/r = -0.92, Oc/r = 1.19. The measured 
curve was multiplied by a factor of 5150 to make the ampli- 
tudes equal, (b) Scattering rate as a function of cooling laser 
detuning for typical experimental conditions: fic/F = 1.19, 
J7r/r = 1.8, and repumper blue-detuned. 



difficult to determine exactly from such a half spectrum 
with possible changes in ion fluorescence near line center, 
but can be seen to be < 2r. 



IX. CONCLUSIONS 



In this paper, we have studied the different types of 
dark states that can occur in a single trapped ^^Sr''' ion 
due to the cooling and repumping laser fields. Efficient 
Doppler cooling and fluorescence detection of the ion re- 
quire a high scattering rate, which is only achieved if dark 
states are prevented from forming. Dark states within the 
55^51/2 ground state can be prevented simply by using 
a linearly polarized cooling beam. Dark states between 
5s "^Sxji ground state and 4d ^1)3/2 metastable state sub- 
levels can be tuned away by using a positive repumper de- 
tuning (assuming that the cooling laser detuning is neg- 
ative for Doppler cooling) when the lasers have narrow 
lines and modest intensities so that the CPT resonances 



are not heavily broadened. On the other hand, the cross- 
correlation between the two lasers and their linewidths 
affect the dephasing rate of the ground-metastable co- 
herences. This means that if the two lasers are uncor- 
related and the sum of their linewidths is of the order 
of the natural linewidth of the ion, the contrast of the 
ground-metastable dark states will vanish. Dark states 
within the metastable state are thus the main problem. 
These can be destabilized using a magnetic field or, in 
applications where strong fields cannot be tolerated such 
as ion clocks, by modulating the repumper polarization. 

Our results regarding the SP linewidth and optimum 
parameters essentially agree with those of Berkeland and 
Boshier [17]. Power broadening by the cooling laser be- 
comes significant for Jlc/r > l/-\/3 « 0.6. However, for 
many practical applications it might be desirable to in- 
crease the Rabi frequency to f2c/r ~ 1, which increases 
the scattering rate by more than a factor of 2, while the 
linewidth increases only by approximately 15%. If the 
repumper Rabi frequency is lower than the cooling laser 
Rabi frequency, there is significant broadening due to 
optical pumping into the metastable state. On the other 
hand, both lasers cause power broadening of the CPT 
resonances. A too high fir value causes the CPT reso- 
nance to extend to the red side of the line. Depending on 
the exact values of J7c and a repumper Rabi frequency 
of Or « (1 . . . 1.5)0c gives the maximum scattering rate, 
while the CPT resonances remain sufficiently narrow that 
they can be tuned away by a positive repumper detuning 

8, « r/2. 

We have also analyzed the effect of ion motion for cases 
where its amplitude is low, i.e., the ion is sufficiently 
cooled and does not undergo excess micromotion. We 
found that at a 5 mK temperature, high-contrast CPT 
resonances, associated with a low ground-metastable co- 
herence dephasing rate 7g,m, i-e., very narrow linewidth 
or cross-correlated lasers, are rounded off by the ion mo- 
tion. For higher dephasing rates 7g.m ^ 0.04, at lower 
temperatures, or further from CPT resonance, the effect 
of ion motion can be neglected. 

The theoretical results have been compared to exper- 
imental data from the ^^Sr"*" endcap trap at the Na- 
tional Research Council of Canada and the agreement 
was found to be very good considering the uncertainties 
of the experimental parameters. The results confirm that 
the system can be well described by the current theoret- 
ical framework and can shed light on the full optimiza- 
tion of single ion fluorescence and laser cooling for sin- 
gle trapped ions applied to atomic frequency standards 
or test systems for quantum manipulation studies. The 
rich variety of observed behaviors can now be understood 
and optimal parameters chosen depending on the desired 
operating conditions needed. It is anticipated that these 
results will certainly aid in providing the optimal high 
fluorescence rates needed for efficient detection of single 
ions and the low kinetic temperatures from laser cooling 
needed for a perturbation free nearly isolated quantum 
system. 
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